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O^' equations of motion in an explicit form. We determine the sectors of solutions that 
{^j;), correspond to different spacetime geometries, and prove that the parameters employed 



We carry out the quantization of the full type I and II Bianchi models following the non- 
perturbative canonical quantization program. These homogeneous minisuperspaces 
are completely soluble, i.e., it is possible to obtain the general solution to their classical 



to describe the different physical solutions define a good set of coordinates in the phase 
space of these models. Performing a transformation from the Ashtekar variables to 
this set of phase space coordinates, we endow the reduced phase space of each of 
these systems with a symplectic structure. The symplectic forms obtained for the 
type I and II Bianchi models are then identified as those of the cotangent bundles over 
£^ ^-j X S"^ X 5"^ (modulo some identification of points) and £^ _|_^ x S^, respectively, 
with £^ the positive quadrant of the future light-cone. We construct a closed 
*— algebra of Dirac observables in each of these reduced phase spaces, and complete the 
quantization program by finding unitary irreducible representations of these algebras. 
The real Dirac observables are represented in this way by self-adjoint operators, and 
the spaces of quantum physical states are provided with a Hilbert structure. 
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I. Introduction 

Quantization of gravity is one of the main obstacles that Modern Physics has to face in 
order to obtain an unified treatment of all physical interactions. To achieve this goal, 
a new and promising approach, the non-perturbative canonical quantization program, 
has been developed systematically over the last years [1-3] . Even though this approach 
has succeeded in solving a variety of physical problems [3] , the implementation of this 
quantization program for the full theory of gravity remains incomplete. 

To apply the canonical quantization program to a given theory, one first selects an 
over-complete set of complex classical variables in the phase space of the system that 
is closed under the Poisson-brackets structure [4] . This set is promoted to an abstract 
*— algebra of quantum elementary operators, with the complex conjugation relations 
between classical variables translated into *— relations in this algebra. The next step 
consists of choosing a complex vector space, and finding on it a representation of the 
*— algebra of elementary operators. The kernel of all the operators that represent the 
constraints of the system provides us with the space of quantum solutions to all physical 
constraints in the considered representation, i.e., with the space of quantum states. 
At this stage of the quantization program, one should determine a set of "real" Dirac 
observables for the theory, that is, a set of operators which correspond to real classical 
variables and such that they commute weakly with all the quantum constraints [4]. 
If this set is suflBciently large, one can fix uniquely the inner product in the quantum 
physical space by requiring that the "real" Dirac observables are promoted to self- 
adjoint operators (reality conditions) [2,3,5]. The space of quantum states is endowed 
in this way with a Hilbert structure. Finally, to extract predictions from the quantum 
theory so constructed, one has to supply the obtained mathematical framework with 
a physical interpretation. 

Owing to the great complexity of the general theory of gravitation, both the com- 
plete space of physical states and a satisfactory set of Dirac observables for gravity 
are still to be determined in the new variables formalism. The analysis and quantiza- 
tion of minisuperspace gravitational models, following the non-perturbative canonical 
program, can be at this point a helpful way to develop some insight into the kind of 
problems, methods, and techniques that are involved in the quantization of the full 
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theory of gravity. In addition to the lessons one can learn by applying the quantiza- 
tion program to simple models, it is clear that obtaining consistent quantum theories 
for minisuperspace models of cosmological interest is physically relevant by itself, as 
it enables us to address cosmological problems quantum mechanically. It is therefore 
not surprising that the recent literature contains a considerable number of works on 
canonical quantization of minisuperspace models [6-11]. 

A special attention has been paid, in particular, to the quantization of the Bianchi 
models [7-10]. These are spatially homogeneous spacetimes which admit an isometry 
group that acts transitively on each leaf of the homogeneous foliation [12,13]. Never- 
theless, the analysis on Bianchi models, has been restricted almost entirely to the case 
of diagonal models [9-11], i.e., models in which the metric is purely diagonal. Even 
if this reduction is completely consistent [11], it would be desirable to carry out the 
whole quantization program for the full non-diagonal Bianchi models. In this way, 
one could also study the role played by the extra non-diagonal degrees of freedom in 
the quantum version of these systems, and discuss the implications of the diagonal 
reduction from the quantum point of view. 

In this paper, we will analyze in detail the full non-diagonal type I and II Bianchi 
models. These systems are completely soluble, that is, one can obtain the explicit 
expressions of the general classical solution for both of these models. We will prove 
that the parameters that appear in the general solution define a good coordinatization 
of the phase space in these two models. The use of the non-diagonal degrees of freedom 
turns out to be decisive in determining the ranges of the introduced phase space 
coordinates. In fact, the reduction to the corresponding diagonal models would lead 
us to different conclusions about the range of the coordinates that describe the diagonal 
degrees of freedom. This is essentially due to the fact that, in the diagonal case, one 
can consistently consider the spatial directions in the homogeneous foliation as fixed 
once and forever. In the non-diagonal case, however, the requirement of analyticity in 
the introduced coordinatization of the phase space obliges us to deal exclusively with 
different homogeneous three-geometries, which are invariant under any interchange of 
the spatial directions that are not preferred by the symmetries of the model. 

On the other hand, using the explicit expressions of the classical solutions, it is 
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possible to endow the reduced phase space of each of these models with an analytic 
symplectic structure. The symplectic forms obtained in this way can be interpreted 
as those associated to real cotangent bundles over some specific reduced configuration 
spaces. Following the canonical quantization program, wc construct, on each of these 
reduced phase spaces, an over-complete set of classical variables that commute weakly 
with the constraints of the system and form a closed Lie algebra with respect to the 
Poisson-brackets structure. This Lie algebra can be identified as the algebra of the 
Dirac observables of the theory. We will then choose a vector space and find on it 
an unitary irreducible representation of the algebra of observables, so that the real 
Dirac observables are represented by self-adjoint operators [14]. The Hilbert spaces 
determined by this procedure provide us with the spaces of quantum states of the 
studied models. 

The outline of this paper is as follows. In Sec. II we introduce the class of mini- 
superspace models on which we will concentrate in this work, and the main formulas 
needed to carry out our analysis in the new variables formalism. The general classical 
solutions for type I and II Bianchi models are obtained in Sec. III. In Sec. IV we study 
the symplectic structure of the space of physical solutions for Bianchi type II. The 
symplectic form in this space is written in terms of the different parameters contained 
in the general solution. We have to prove then that the chosen set of parameters 
defines a good coordinatization of the phase space of the model. This is the subject of 
Sec. V. Our analysis for Bianchi II is generalized to the non-diagonal type I Bianchi 
model in Sec. VI. Sec. VII deals with the non-perturbative canonical quantization of 
these two models. In Sec. VIII we present a different approach to the quantization of 
Bianchi type I, using the symmetries that are present in this model at the classical level 
to find a complete set of generalized "plane waves" which span the space of quantum 
states. We also show that the two quantum theories constructed for type I arc unitarily 
equivalent. Finally, we summarize the results in Sec. IX, where we also include some 
further discussions. 

II. Bianchi Models 

Bianchi models are spatially homogeneous spacetimes [i.e., they can be foliated by 
three dimensional Riemannian manifolds) which admit a three dimensional isometry 
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Lie group G that acts simply transitively on each leaf E of the homogeneous foliation 
[12,13]. As a consequence, there exists for each of these models a set of three left- 
invariant vector fields Lj on E which form the Lie algebra of the group G: 

[Li,Lj] = C''ijLk , (2.1) 

where jk are the structure constants of the Lie group. Dual to the the vector 
fields Lj, one can introduce a set of three left-invariant one-forms which satisfy the 
Maurer-Cartan equations 

dx' + lc'jK X'^X"" =0 ■ (2.2) 

If the trace C^/j of the structure constants is equal to zero, the Bianchi model 
is said to belong to Bianchi class A. For this class of models, the spacetime admits 
foliations by compact slices [15]. We will restrict ourselves to this case hereafter. 

The structure constants for the class A Bianchi models can always be written in 
the form [12,16] 

JK = ^jklS^^ 1 (2.3) 

with ejKL the anti-symmetric symbol, and S^^ a symmetric tensor. Further classifi- 
cation of the class A Bianchi models is defined with respect to the signature of the 
symmetric tensor S^^ [12,16]. The type I and II Bianchi models, the only ones we 
will consider in this work, are characterized by the signatures (0,0,0) and (0,0,-1-), 
respectively. Thus, the structure constants for Bianchi I are given by C^j^ = 0, while 
for Bianchi II jk = ^sCsjk- 

In the new variables formalism, one starts by introducing the triads e^ai.^) on a 
three-manifold S, where z = 1, 2, 3 is a spatial index and a = (1), (2), (3) is an S0{3) 
vector index which is raised and lowered with the metric rjab = (1)1)1)- The inverse 
metric on E can be written in terms of the triads as g^^ = e*ae^°'- The Ashtekar 
variables {S'^a, ^j^) are defined as the densitized triad and the spin connection [3]: 

£\ = idetg)h\ , A" = r,'^(e) -iK,- , (2.4) 

with Fi" the SO{3) connection compatible with e^a, and Ki"" the triadic form of the 
extrinsic curvature. 
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For the class A Bianchi models and S a compact manifold, one can always perform 
the following transformation of variables [17]: 

{Ai''{x,t),£\{x,t)) ^ {Aj''{t),At{x,t),E\{t),Si{x,t)) (2.5) 

where t is the time introduce by the homogeneous foliation, a; is a set of coordinates 
on the leaf S, and 

= ^ / d'x \x\ Ai^Lr' , At = Ai^ - Ai^x'i , (2-6) 
E\ = U d^xS\x'i, S: = S\-E\Li^\x\. (2.7) 

In Eqs. (2.6,7), = X^idx^^ Lj = Lj^di, \x\ is the determinant of x^ i cind 

= J.^d^x\x\- The reduction of the dynamical degrees of freedom in the class A 
Bianchi models is accomplished by imposing Af = = 0. In this way, one is left only 
with a finite number of degrees of freedom, given by {Ai''{t),E^ait)). The variables 
(Ai^-it), E\{t)) form a canonical set [11,17]: 

{Ai''{t\E\{t)} = i5j5t . (2.8) 

In terms of them, the constraints for the class A Bianchi models adopt the following 
expressions [17]: 

Ga = eab'^Ai'^E', , (2.9) 
Vi = C'ijAk'^E^ , (2.10) 
S = ea'^i-C'^ijAK'' + e'^deA/Aj^) E\e\ , (2.11) 

where Qai Vj, and S denote, respectively, the Gauss law, the vector constraints and 
the scalar constraint. 

On the other hand, using the left- invariant one-forms x^ i the spatial three- metric 
on E can be written as ^ = QuX^X^ i with E^ a related to the inverse of gu by means 
of 

g'\deiE) = E\EJ'' , (2.12) 
and det£^ denoting the determinant of E^ a- 
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Finally, the connection Aj"' admits the decomposition 



Ai"" = Ti''{E)-iKi'' , (2.13) 

where the SO (3) connection F/" and the triadic extrinsic curvature Kj°' can be deter- 
mined through the formulas 

ri- = e\''{^C''ji{E-y^E\-^C''ji{E-y^{E-yjE\E\) , (2.14) 
Ki" = KijE\{detE)-^ , (2.15) 

with {E-^^j the inverse of E^ a and Kjj the extrinsic curvature associated to the 
metric gij. 

III. Bianchi Types I and II: Classical Solutions 

In the rest of this work we will restrict our attention to the type I and II Bianchi 
models [12,18]. It is well known that, for these models, one can always reduce the 
geometrodynamic initial value problem to the diagonal case [19,20]. We will now 
briefly review the argument that leads to this conclusion. 

Suppose that, for either of these two models, we begin by considering a certain set 
of left-invariant one- forms i for which the three- metric is given by gij. Any other 
set of left-invariant forms will be related to x^ by a transformation x^ = jx"^ 
that maintains the symmetries of the model. The structure constants must, therefore, 
remain unchanged under the transformation defined by M^j, 

C'jK = {M-y ^C"^ pqU"" jM^ K . (3.1) 

The inverse (M~^) must always exist, since is a set of three linearly independent 
one-forms. 

For Bianchi type I, jk = 0, and condition (3.1) is empty. In this case, any 
invertible matrix M e GL(3,IR) defines a permissible transformation. Thus, for any 

geometrodynamic initial value data [g^jj^K^u) in the set x^ (with g^jj a positive 
definite metric), we can perform a transformation with a matrix M G '^qq^^^ such 
that, in the new set of one- forms, the initial metric takes the value g^j — Sjj. Then, 
using a transformation under SO{3), we can bring the initial extrinsic curvature K'^ij 
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to the diagonal form, without altering the identity value for the initial metric [19,20]. 
Since the diagonal ansatz is compatible with the dynamics of the type I Bianchi model, 
we conclude as a corollary that any geometrodynamic classical solution for type I can 
be expressed as 

9IJ = {M')i''gDM''j , (3.2) 

where the metric (725 is a classical solution for the diagonal case, M^j is a constant 
invertible matrix, and (M*) denotes the transpose of M. 

Let us consider now the type II Bianchi model. For this model, jk = S^^zjk-, 
and condition (3.1) implies that M^j must be of the form 

M = f ^ 1 ' (3-3) 

V Mf M| detM / 

with M. e GL(2,IR). Given any initial value data {g^jj,K^ij), defined in the set 
of one-forms j we can always carry out a transformation under a matrix of the 
type (3.3), with M. G ^§o(^, such that, in the new set = M^jx"^, the ini- 
tial metric is equal to g'jj = Sjj + (g^ — l)(^/3(^j3, and the extrinsic curvature is 

Using expressions (2.12-15), the Gauss law constraints (2.9) and the two non- 
empty vector constraints in (2.10) for Bianchi type II, it is possible to see that, for 
positive definite metrics g^j of the form that we have obtained, K^^ and K^2 must 
vanish if (^/j, Kjj) is an admissible set of initial value data. With these conditions, 
it is clear that Kjj can be brought to diagonal form by a transformation of the type 
(3.3) with M G S'0(2) and M^i = M^2 = [19]. Under such a transformation, the 
initial metric g^j remains unchanged. Since the diagonal case is consistent with the 
dynamics of the model, we conclude, as for Bianchi I, that any Bianchi type II classical 
solution can be written in the form (3.2), with M an invertible constant matrix of the 
type (3.3). 

Therefore, to get the general solution in geometro dynamics for the Bianchi types I 
and II, it suffices to find the classical solutions for the corresponding diagonal cases. 
These solutions can in fact be obtained from the analysis of the diagonal Bianchi 



models made by Ashtekar, Tate and Uggla in Ref. [9]. Parallelling their notation, we 
introduce the following parametrization for the diagonal metric g^, D = 1, 2, 3: 

In the Misner's gauge [21], defined by the lapse function = 12{detgD)^ 
= 12 6^*^^^°"^+-*, the dynamical equations for Bianchi type I adopt the simple ex- 
pression 

Po = -Po , P+^P+, P-^P-, (3.5) 

where the dot denotes time derivative and (po,P+,P-) are three real constants which 
are related through the scalar constraint [9] 

S <x -po^ + +pJ = . (3.6) 

Integrating the equations of motion, we arrive at a diagonal metric and a lapse 
function ^ 

g^ = e-2Vsipo+P+-P-)t ^ g^ ^ ^-2V3ipo+P++P-)t ^ g^^^2V3p+t^ (3_7) 

N = 12e-^^'^P°+P+^^ . (3.8) 
Finally, from (3.6), we can restrict our analysis to non- negative po given by 

Po = Vp+'+P-' • (3.9) 

The solutions corresponding to negative po can be obtained from those with po > by 
changing the sign of the time parameter t that defines the evolution and hipping the 
signs of p+ and p-. All the different physical solutions {i.e., solutions with different 
spacetime geometries) are contained in the sector po > 0,p^,p-,t G IR [9]. We will 
thus restrict ourselves to this range of the parameters appearing in (3.7-9). 

To find the classical solution for Bianchi type II, one needs to perform a canonical 
transformation that mixes with its canonical momentum [9]: 



^ The integration constants for Po, P+ and /3_ can be absorbed by a translation of 
the origin of time and a redefinition of the matrix M that appears in (3.2). 
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with defined as a strictly positive variable. In Misner's gauge, the dynamical equa- 
tions (3.5) and the scalar constraint (3.6) are still valid for the type II diagonal model, 
with the substitution (/3_|_,p_|_) (/9_|_,p_|_), and (po,P+,P-) three real constants. It 
is then straightforward to derive the general expression for the diagonal metric in the 
classical solutions 

^ ,-2VsiPo-P-)t 2V3cosh(2v/3p-+t) 

= e-2V3(po+P-)t 2v/3cosh(2V3p+t) ^ ^^^^^^ 

P+ 

~ 2V3cosh(2V3p+t) ' 
where t is the time coordinate defined by the lapse function 

N = i2e-'V3pot /2V3cosh(2v/3p+t)x h 
\ p+ J 

and 



p+>0, p-eTR, po = a/p^^^^P^ > . (3.13) 

Using the explicit expressions (3.7,8) and (3.11,12) for the metric and the lapse 
function in the diagonal type I and II Bianchi models, the relations (2.12-15) and (3.2), 
and the formula 

Kij = ^gij , (3.14) 

valid in the gauge in which the shift functions vanish [13], one can easily compute 
the general form of the Bianchi I and II triads and spin connections in the physical 
solutions, restricted to the sector that corresponds to positive definite metrics. The 
result can be written in the compact notation 

E\ = detM{M-^y ^EdTZ"^ , (3.15) 
yl/« = (M^)/^ ^ 7^"■^ , (3.16) 

where 7^"^ is a general complex orthogonal matrix, and the sum over D = 1, 2, 3 is 
implicitly assumed. 
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For Bianchi type I, the 



matrix M belongs to GL(3,IR) and 



(3.17) 



{Po + P+ 



P-), ^1 



{PO+P++P-), 



— 1 



(3.18) 



Ui = 



4^3 



4^3 



4V3- 



with Po given by (3.9). 

In the case of Bianchi II, the matrix M appearing in (3.15,16) must be of the form 
(3.3), and {Ed,ojd) can be expressed as 



where {po,p+,p-) must satisfy the restrictions (3.13). 

Even if the classical solutions that we have found correspond to positive definite 
metrics, it is clear from (3.15) that one can always reach degenerate metrics in the 
limits in which either Ei,E2,E3 or detM vanish. It is only in this sense that the 
degenerate solutions are included in our analysis. 

IV. Bianchi Type II: Symplectic Structure of the Space of Solutions 

In this section, we will concentrate our attention on the type II Bianchi model, studying 
the structure of the space of physical solutions. 

Let be the set of left-invariant one-forms for which the metric gu is diagonal. 
For Bianchi type II there always exists a preferred one-form selected by the sym- 
metries of the model, since the structure constants are given by jk — S^^eljk, 
with 5"^"^ the only non- vanishing component of 5''^'^. However, the spacetime geome- 
tries remain obviously unaltered under the interchange of and x^. Two classical 
solutions which are related by the interchange of indices 1=1 and 1 = 2 should then 
be identified as the same physical solution. ^ From the expressions (3.19,20), this 

^ Nevertheless, one can neglect this identification by considering x^ and x^ as two 
preferred one-forms. The space of solutions will then correspond to a different theory, 
in which the discussed symmetry is not present [9]. 



E, = e--^(P°+^-)S E, = e-^/^(Po-P-)t ^ ^3 = 2v/3cosh(2v/3p+t) 



P+ 
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interchange of indices can be realized as a flip of sign in p_ . Therefore, we can restrict 
ourselves only to non-negative parameters p_ > 0, so that each physical solution is 
considered only once. 

There is still some redundancy left in the classical solutions of our model. This 
redundancy comes from the fact that, if ^ is a matrix that satisfies 

A'gDA = gD (4.1) 

for all diagonal metrics go, and such that AM is of the form (3.3) for every matrix M 
of that form, the classical metrics (3.2) associated to the collection of matrices AM 
turn out to be identical [22]. The conditions imposed on A define a discrete group of 
four elements 

{A, = (1, 1, 1),A2 = (1, -1, -1) , ^3 = (-1, -1, 1),A^ = (-1, 1, -1)} , (4.2) 

where (a, 6, c) denotes the ordered set of diagonal elements of A, and all the non- 
diagonal elements are equal to zero. Using the invariance of the physical solutions 
under multiplication of M by we can choose a positive determinant for the matrix 
M. appearing in (3.3). We will thus restrict in the following to the case detAl > 0. 
Note that, however, we have still to identify the classical solutions corresponding to 
M and A3M, since multiplication by ^3 conserves the sign of detM.. We will return 
to this point later in this section. 

Once we have determined the physically different classical solutions, we proceed 
to show that the space of solutions is endowed with a symplectic structure. We begin 
with the symplectic structure in the Ashtekar formalism: 

in = dAi''AdE\ . (4.3) 

Substituting Eqs. (3.15,16) in this formula, we arrive at the expression 

in = d{uidetM) A + d{detM) - {dMM'^ 

- uidetM {dMM-y J A{dMM-^Y J , (4.4) 

with {dMM-^y J = dM^ q{M-^)^ J and d{detM) = detM Tr(dMM-i). The complex 
orthogonal matrix TZ"'^ that is present in Eqs. (3.15,16) disappears completely from 
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the symplectic form (4.4), as it represents only the gauge degrees of freedom associated 
to the Gauss law constraints (2.9). On the other hand, taking into account that the 
matrix M is of the form (3.3), with detAd > 0, we can always decompose M in the 
following product of matrices 

M = MdMsMtR 

/aOO\/l 0\/lOO\/cos^ sin^O\ 
= 06 1 o LzlO -sin^ cos^ , (4.5) 

\0 ab J \Mf M| 1/ \0 1/ \ 1/ 

with a, 6 > 0, Mf , M| > z eJR, and 9 E S^. The diagonal matrix Md can be absorbed 
into the diagonal part of the classical solutions, (£■£>, cud), by means of the redefinitions 

El = ah^Ei , E2 = a^bEi , E3 = abE^ , Qd = {abfuo ■ (4.6) 

Introducing then the notation M = M^MtR, Eq. (4.4) can be rewritten as 

in = d{uji) A (d{\iiEi) - {dMM-y - o)i{dMM-y j A {dMM'^y j , (4.7) 

where we have employed that {dMoM-^DY j = S^jiidMoM'^ dY j- 

Let us define S = MtR, so that M = M3S. It is straightforward to compute that 

idMM-'Yj = diM^S] + d{Ml)5i5] + {M:,Y QidSS'^)"^ p{M^Y j ■ (4-8) 

Then, using that {dSS'^Ys = and 

{M^YiMMsYj = Di + Di6i{M^6'Q + mIS^) (4.9) 

for D any diagonal matrix, we conclude that Eq. (4.7) is still valid with the substitution 
of S for M. In this way, the matrix (M3) drops from the symplectic form ft, implying 
that Mf and Mf do not correspond to real degrees of freedom. 

In the parametrization chosen in (4.5), {dSS~^) takes the explicit expression 

/ -zd9 de 0\ 

{dss-^) = dz-{z^ + i)de zde , (4.10) 

\ 0/ 
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from which it follows that 



in = d{o)D) A dilnEo) + d{{ai - Sa)^) A d9 . (4.11) 

Substituting Eqs. (4.6) in this formula, with (EdjLOd) defined by means of (3.19,20) 
and (3.13), we arrive after some calculations at a symplectic structure of the form 

n = dU+ AdX + dn_ AdY + dOAdZ , (4.12) 

where 

n+ = {aifp+, n_ = (o(,)V-. (4.13) 

Expression (4.12) provides us with the symplectic form for the space of solutions to all 
physical constraints, i.e., the reduced phase space. Nevertheless, in obtaining this sym- 
plectic structure, we have implicitly assumed that the parameters (p+,p_,po) t-, o-i b, z, 
9, Mf, Mf ) good coordinates in the space of physical solutions to all but the scalar 
constraint; that is, that the transformation from the triad and the spin connection to 
the given set of parameters is analytic in the sector of the phase space covered by these 
solutions. We will prove that this is indeed the case in Sec. V. 

On the other hand, introducing the notation Hq = {ab)'^po, the scalar constraint 
(3.6) implies that (n+,n_) are a set of coordinates on the light-cone: 

-Ho^ + + UJ = . (4.15) 

We point out that, in our model, p^, a and b are positive quantities, and that we have 
to restrict ourselves to positive po and non-negative p- in order to deal exclusively 
with different physical solutions. With these restrictions, {X,Y,Z) defined in (4.14) 
run still over the whole real axis, for z eJR. We can then interpret the two first terms 
on the right hand side of (4.12) as the symplectic form of the cotangent bundle over 
'^(+ _|_) 1 the positive quadrant of the future light-cone. 

In the last term of (4.12), the angle 9 belongs to S^. We recall, however, that, 
from our previous discussion, there is still some redundancy to be removed if we want 
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to consider only physically different classical solutions. This redundancy corresponds 
to the identification of the matrix M appearing in (4.5) with that obtained by multi- 
plication on the left by As = (—1, —1, 1). As a consequence, the solutions associated 
with the parameters {9, Mf , Mf) and {9 + tv, — Afi , —M2) are physically identical, and 
we may restrict our analysis to the interval 9 G [0,7r). Note, nevertheless, that Mf 
and Ml are not real degrees of freedom (we can always go to the gauge in which 
Mf = Ml = 0)) that the coordinates (n+, n_, X, F, Z) remain unaltered under 
the transformation (6*, Mf, M|) {9+7r, -Mf , -M|)- Therefore, the identification of 
physical solutions under that transformation obliges us to identify also the boundaries 
^ = and ^ = TT of the reduced phase space. In this way, the angle = 2^ turns out 
belong to S^, and the term d9 A dZ = dO A d{^) in (4.12) can be interpreted as the 
symplectic form of the cotangent bundle over . We thus conclude that the space 
of physical solutions in Bianchi II presents the symplectic structure of the cotangent 
bundle over the reduced configuration space _|_^ x . 

V. Bianchi Type II: Analyticity of the Coordinatization of the Space 
of Solutions 

We want to prove that the parametrization employed to describe the classical solutions 
for Bianchi type II (with po regarded as independent of p_ and p^.) defines an analytic 
coordinatization in the space of physical solutions to all but the scalar constraint, 
in the sense that the transformation from the triad and the spin connection to the 
chosen set of parameters is analytic in the whole region of the phase space covered 
by the different physical solutions that we have considered. We note first that the 
transformation from the triad and the spin connection to the three-metric and the 
extrinsic curvature is analytic in the sector of the Ashtekar variables that corresponds 
to non-degenerate metrics. It will suffice then to show that the matrix M, appearing 
in (3.2), and the parameters {j>j^,p-^pQ,t), contained in qd^ depend analytically on 
((7/j, Kij) in the region defined by detA1,p-|-,p_ > 0, which contains all the different 
physical solutions. 

Let us introduce a matrix M of positive determinant (and thus invertible) such 
that it satisfies the conditions 

{{M-'Y) gPQ{M-')^ J = 5jj + {gs - l)^!^} , (5.1) 
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(5.2) 



It is clear from our previous analysis of the type II Bianchi model that one solution to 
Eqs. (5.1,2) is provided by 

M= [igi)K{g2)Ki9i92)^)M , (5.3) 

where we have used a similar notation to that displayed in (4.2), M is the matrix 
appearing in (3.2,3), and (^1,52) are given by (3.11) in our parametrization. From 
Eqs. (3.11,12) and (3.14), one can compute also the explicit expressions of and A/ 
for the solution (5.3): 



gs 



4A^( f+ -)', Xi=X{-po+P-+p+ta,nh{2V3p+t)), (5.4a) 

^cosh(2v3p+t)^ 

A2 = X{-po -P-+P+ tanh(2v^p+t)) , A3 = -Xg3P+ tanh(2v^p+t) , (5.46) 

X=—^{^ ) • (5-4c) 

4^3 V2^cosh(2v^p+t)>' 

Since g\ and g^-, given by (3.11), are strictly positive, Eq. (5.3) defines M analyti- 
cally in terms of (p+,p_,po5 ^5 -^)- All we have to prove then is that (p+,p_,po5 ^5 -^) 
can be obtained analytically from (^f/j, Kij) in the region of the phase space that we 
are considering. On the other hand, it is easy to check that the matrix M in (5.3) 
is of the form (3.3), provided that M is of this form. Moreover, adopting a parallel 
notation to that introduced in Eq. (3.3), one can see that detTW > if detA^ > 0. 
We will thus concentrate in the rest of this section on matrices M of the type (3.3), 
for which detAi > 0, and such that they verify Eqs. (5.1,2). 

It is convenient to employ the following decomposition for the matrix M: 

M = RMdMtMs , (5.5) 

with R, Mdi Mt, and M3 parametrized as their respective counterparts in Eq. (4.5) 
with the substitutions {a,b, Mf, M2, z, 9) — > {d,b, Mf, M2, z,6), so that a, 6 > 0, 
z, Mf, M| G IR and 9 E S^. Let us show then that all the elements of these matrices 
present, through Eqs. (5.1,2), an analytic dependence on {gij,Kij) in the region of 
physical interest. 
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We first determine M3 by the conditions 

^13 = ^23 = , with gij = {{M^y)^''gpQ{M^'fj , (5.6) 

from which one obtains that Mf = gis/gss and = 912/933- M3 is thus analytic in 
gij, since gu is a positive definite metric, and therefore ^33 > 0. Let us define then 
the two-by-two matrix h constructed with the two first rows and columns of g, and, si- 
milarly, the matrix K, obtained from the extrinsic curvature K = (M^^yK(M^^). ^ 
We can use the degrees of freedom a, b, and ^ in M to diagonalize h to the identity 

(M-r^yiM^y h iM^'XMr^') = I , (5.7) 

where M.d and M.t are the two- by- two matrices constructed from Md and My by 
the procedure explained above. Eq. (5.7) fixes a, b, and z uniquely, for a and b strictly 
positive, 

a= ^ — , 6=(/i22)2, z=^. (5.8) 

^ /i22 ^ "22 

From the expression (5.8) we conclude that a, b, and z are analytic in h (and then in 
g) for h positive definite. 

We can now diagonalize £ = (A^~^)*(A^^^)*^(A4^^)(>(^^) by the 50(2) trans- 
formation given by the matrix TZ contained in R. This transformation leaves invariant 
the value h = I reached in (5.7). We are thus left with the eigenvalue problem 

Fij = 1Zj'')CpQ{1Z'fj - X161J = . (5.9) 

The element F12 of the system of equations (5.9) defines 6 as an implicit function of K,. 
Taking into account that, from our previous analysis, K, is analytic in (gu, Ku), all 
we have to prove is that the implicit dependence of 6 on K., imposed by Eq. (5.9), is 
analytic. Let us suppose that {9''\ IC^j) is a particular solution to the equation F12 = 
that determines as a function of IC. Then, 6{IC^) = 9^ defines an analytic germ 
[23] around IC^ which can be continued analytically as long as d0Fi2(^9{)C), K.) ^ 



^ From our discussion of the type II Bianchi model in Sec. Ill, it follows that, given 
the form of g, K must satisfy K12, = K23 = if {gij, Kij) corresponds to a classical 
solution to all but the Hamiltonian constraint. 
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[23,24] (in particular, if $^ and iC'jj are real and 0{]C) can be continued analytically 
around ^/j, it is possible to show that d{K,) remains real for }C eJR). After a simple 
computation, one arrives at the identity 

deF^2 = X2-Xi , (5.10) 

so that the obtained solution 9{K,) depends analytically on K, as far as Ai 7^ A2. ^ 
Note that, in our parametrization, Ai = A2 only if p_ = (see Eqs. (5.4)). Since 
we are restricting our attention to the sector of positive definite metrics and extrinsic 

curvatures for which p- > 0, we conclude that 0{IC) is analytic in that region, except 
at the boundary po = 0. ^ With this caveat, the matrix M turns out then to be 
analytic in the sector of metrics and extrinsic curvatures associated to the physical 
solutions analyzed in Sec. IV. 

Employing this result, Eqs. (5.1,2) provide us with and A/ (/ = 1,2,3) as 
analytic functions of {gij, Kjj), for (M~^) is always well-defined, and is analytic if M 
is analytic. We can use now the explicit expressions (5.4) to determine (p_|_, p_, po, t) as 
analytic functions of ^3 and A/, and, therefore, of {gu, Ku). The relations obtained 
are 

P-=P+^, Po=P+^, (5.II0) 

-^3 -^3 

= 2V3Fs[±e-'^'""'-'^^f , (5.11c) 



where 



= J(Ai - A2) , F2 = - J(Ai + A2) , (5.12a) 

2 6-3 2 



^ Eq. (5.9) can be interpreted as an eigenvalue problem. Each row of the matrix TZ 
can be identified as an unit eigenvector of the matrix IC. When Ai = A2, the eigenvalue 
problem is degenerate: any unit vector is an eigenvector of )C. As a consequence, TZ is 
ill-defined in that case. 

^ It is clear then that, without a restriction of this type, p- would have not been a 
good coordinate in the phase space of the full type II Bianchi model. 
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For real positive definite metrics (so that > 0) and real extrinsic curvatures, all the 
functions F„ (with n = 1, 4) are analytic in (^3, A/), and F3 and F4 turn out to be 
positive. Then, {p^,p-,pQ,t), given by (5.11), result to be analytic in {gu, Kjj). 

From our previous discussion, it follows that the matrix M, defined by means of 
(5.3), is analytic in {gij,Kij) in the sector of solutions to all but the Hamiltonian 
constraint that we are studying. The only point that remains to be proved is that the 
specific parametrization (4.5), employed for M in Sec. IV, is analytic with respect to 
its dependence on the elements of M. Using this parametrization, and the notation in 
(3.3), we have that 

KA_(M\ M^2\ _ ( a OW cos^ sin^\ . . 

•^"l^M^i M\)-\hz h)\-sin9 cose)' ^^■^'^> 

Ml = ah {mI cos e + M| {z cos ^ - sin ^)) , (5. 14a) 

Mf = ah^Ml sin e + Ml{z sin e + cos 9)^ , (5.146) 

where a and h are strictly positive, and det > by construction. From Eq. (5.13), 
one can easily arrive at the explicit expressions of a, 6, and z in terms oi M.M.^. These 
expressions are analytic in the elements of M. for detAl > 0, and allow us to define a 
and h as strictly positive functions: a{M), h{M) > 0. Eq. (5.13) provides us also with 
the relations cos^ = Ml/a{M) and sin^ = M2/a{M), from which we conclude that 
cos^ and sin^ are analytic in Ai, as a{Ai) > 0. ^ Substituting then a,b, z,cos9 and 
sin^ as functions of Ai, the equations in (5.14) determine Mf and Mf as two implicit 
functions of the elements of the matrix M : (Mf (M), Ml (-^))- The Jacobian of these 
equations with respect to Mf and can be computed straightforwardly to be equal 
to a'^{A4)b'^{A4). Therefore, given a particular analytic germ for (Mf,M|), we can 
always continue it analytically to the whole range of matrices M with detAl > (so 
that a{M) and b{M) are strictly positive) [23,24]. 

In conclusion, we have proved that the parameters {p+,p-,po, t, a, b, z, 9, Mf, Ml) 
are good coordinates in the space of non-degenerate real solutions that satisfy all but 
the Hamiltonian constraint, when restricted to the region 

p+,p-,a,b>0, 9 e , po,t,z,Mf,MieM. (5.15) 

^ Indeed, these relations define the angle 9 analytically, because the Jacobian matrix 
of (cos 9, sin 9) with respect to 9 has always a rank equal to the unity. 
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If we want to go to the reduced phase space, we have to impose in addition the 
constraint (3.6), and restrict the range of po to the positive real axis. As a consequence, 
the symplectic structure obtained in Sec. IV results to be analytic everywhere in the 
cotangent bundle over the reduced configuration space £^ x 5^, with its boundary 
excluded. 

VI. Type I Bianchi Model 

Let us proceed now to generalize the study of Sees. IV and V to the case of the 
type I Bianchi model. We begin by analyzing the sector of classical solutions that 
corresponds to different spacetime geometries. 

In the set of left-invariant one-forms which the metric is diagonal, the 

spacetime geometries are invariant under permutations of all ^ because the structure 
constants vanish identically in this model, so that there is no preferred one-form. 
Therefore, the classical diagonal solutions that are related under any interchange of 
the indices / = 1,2,3 must be identified. From expressions (3.17,18), each of the 
planes Hi = p_ = 0, II2 = po + 2p+ — p_ = and lis = po + 2p+ + = divide 
the space (po,P+,P-) G into two regions which can be interchanged under the 
respective permutations of indices: 1 ^ 2, 1 ^ 3 and 2^3. If we take into account 
also the constraint (3.9), which implies that (po,P+,P-) lie in the future light-cone, 
the requirement of considering only different classical solutions may be implemented 
by the following restrictions in the ranges of our parameters: 

p+>0, P->0 Po = \/p+^ +P-^ > 0, (6.1) 

i.e., {po,p+,p-) e >C^+,+)- 

Parallel to the situation in Bianchi type II, we have to identify also the solution 
(3.15,16) for a given matrix M e GL(3,1R) with all other solutions obtained from 
matrices of the form AM, where A is any orthogonal matrix that commutes with the 
diagonal subgroup of GL(3,IR) (see Eq. (4.1)). We can then fix the determinant of 

M to be strictly positive, because either M or —M has a positive determinant, and 
A = (— 1,— 1,— 1) is orthogonal and commutes with the diagonal subgroup. There 
is, however, some redundancy left, because there exist matrices A which conserve the 
sign of the determinant of M. These matrices are given by the discrete group of 
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four elements defined in (4.2) [22]. We will discuss the corresponding identification of 
classical solutions later in this section. Finally, we point out that, given the explicit 
expression (3.7) for the diagonal metric in Bianchi type I, it is always possible to absorb 
the determinant of M by a redefinition of the origin of the time coordinate t. We will 
thus restrict hereafter to matrices M e 5'L(3,1R). 

Let us introduce the following parametrization for the matrices M e SL{3,JR): 

fa 0\ /l 0\ 
M = MdMtR = h O x 1 f i?) , (6.2) 

where a and h are strictly positive, y and z are real and R e S0{2t). Similar to the 
notation for Bianchi type II, we will call 

Ei = —, E2 = ^, E3 = abEs. (6.3) 
a 

Then, it is straightforward to generalize the discussion presented in Sec. IV to arrive 
at a symplectic structure of the form (4.7) for the space of solutions in Bianchi type I, 
with M = MtR and a)/ = oji (note that now detM = 1). 

The symplectic form (4.7) can be simplified as follows. Decompose first the matrix 
M as M = MtR, and define T^-^ = dR^^R-^^. From the expression of Mr in (6.2), 
one can explicitly check that 

Y,d{MTyp{MT-Yi = (6.4) 
p 

for every / = 1,2,3. On the other hand, the matrix of one- forms T^"^ turns out to 

be antisymmetric, since R^^R^^ = 5^'\ and so dR^^R"^^ = —dR^^R^^. Using this 
fact, it is possible to prove that 

dT" = T^^ A T^-^ . (6.5) 

Conveniently rearranging the different terms in (4.7), we arrive at the conclusion 

\Vt = d(ji ^d\n{EI)-d{LIJT^•^) , (6.6) 
Lij = {Mt~\i\pUp{Mt)^ , (6.7) 
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where the indices /, J are raised and lowered with the metric rju = (1, 1, 1) and [ | ] 
denotes antisymmetrization. The first term on the right hand side of Eq. (6.6) can be 
calculated from Eqs. (6.3), (3.17,18) and (3.9). The result is 

(LjiA d\n{Ei) = -^{dp+ AdX + dp_ A dY) , (6.8) 
4v 3 

with 

X = -(^ + 2) ln(ab) , Y = \n(y)-^ Iniab) . (6.9) 

The one-form Tr(— LT) that appears in (6.6) can be interpreted as the pre-symplectic 
structure in SO{3). From the parametrizations for ui and Mt given by (3.18) 
and (6.2), one can explicitly compute the matrix L/j. Introducing the notation 
Lij = sjjkI^ , we find that 

h = {po + 2p++p_)z , (6.10a) 
l2 = - (2p_{xz -y) + {po + 2p+ +p-)yj , k = -2p-x . (6.106) 

Let us parametrize now the matrices R e SO{3) in terms of the Euler angles {a,P,9) 

i?(«, A ^) = i?(i)(a)i?(3) (/?)%) (^) , (6.11) 

where /?(/)(«) is a rotation of an angle a around the axis defined by the direction /, 
and a, 9 e , (3 G [0,7r]. This parametrization is unique and well-defined for all 
matrices R e SO{3), except at /? = and /? = tt. After a short calculation using Eqs. 
(6.10,11), we conclude that the pre-symplectic structure Tr(— LT) can be written as 

LijT^-^ = -^{Fida + F2d(3 + F^dO) , (6.12) 
4v3 

with 

Fi = Zi , F2 = ^2 sin a + ^3 cos a , (6.13a) 

F:i = —I'iSm asm 13 + I2 cos a sm. 13 + li cos 13 . (6.136) 

Therefore, the symplectic form (6.7) has the expression 

Q = -^{dp+ AdX + dp- AdY + daA dFi + d(3 A dF2 + dO A dFs) . (6.14) 
4v 3 
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For a and b strictly positive and x,y,z real, the variables X,Y, Fi, F2 and F3 run 
over the whole real axis IR. As it stands, the symplectic structure (6.14) might be 
interpreted as that corresponding to the cotangent bundle over £^ ^-^ x 5*0 (3). Nev- 
ertheless, we have still to identify the classical solutions obtained from all matrices 
of the form AM, with M G SL{3,M) and A any matrix in the discrete group (4.2). 
In the parametrizations (6.2) and (6.11), the change of M to A2M can be realized 
as the transformation [a,x,y) (a + tt, —x, —y), which leaves invariant all the vari- 
ables appearing in (6.14) except a (see Eqs. (6.10) and (6.13)). In order to consider 
only different physical solutions, we can restrict ourselves to the range a e [0,7r) [22]. 
The boundaries a = and a = tt must be identified in the reduced phase space, 
since, from our previous discussion, the classical solutions for o; = and a = tt are 
physically identical. In this way, a = 2a belongs to 5"^, and we can interpret the sym- 
plectic form da A dFi +d(3 A dF2 as that of the cotangent bundle over the two-sphere, 
5'^, parametrized by the angles (3 and d. The singularities (3 = and /3 = tt of the 
parametrization (6.11) then correspond to the poles of this two-sphere. 

Let us now study the identification of the matrices M and A3M. The interchange 
of these two matrices is performed by the transformation (a, /3, 9, y, z) ^ {n — a, 
TV — + —y, —z), where, we recall, a G [0, tt). We could then impose the restriction 
9 G [0,7r), so that each physical solution is considered only once [22]. Note that, 
under the above transformation, p^,p_, X,Y and F3 remain invariant, while Fi and 
F2 flip their sign. We should thus identify the points {a, (3,6 = 0, Fi, F2, F^) and 
(tt — a, TT — ^ = TT, — Fi, — F2, F3) at the boundaries ^ = and ^ = tt of the reduced 
phase space. Nevertheless, we will adopt a different approach for the quantization of 
the model, leaving 9 to run over the whole of and imposing restrictions on the space 
of physical states associated to the identification of points 

{a,p,9,F^,F2,F3) and (tt - a, tt - ^ + tt, -Fi, -F2, F3) . (6.15) 

We will return to this issue in the next section. 

Finally, the matrices M and A4M are interchanged through the transformation 
(sin a, /3, 9, x, z) (— sin a, tt — /3, ^ + tt, —x, —z), which leaves invariant p^,p-, X, Y, 
and F3, and reverses the signs of Fi and F2. For a G [0, 7r),sinQ; is non- negative, 
and the only possible redundancy left is at a = 0. However, we notice that the 
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corresponding identification of points {a = 0, P,9, Fi, F2, F^) and {a = 0, tt — 
6 + 7r, —Fi, — F2, F3) has already been taken into account by identifying the points in 
(6.15), and a with tt + a. 

In conclusion, the space of physical solutions in the type I Bianchi model has 
the symplectic structure of the cotangent bundle over the reduced configuration space 
X 5^ X S'^. To consider only different physical solutions we still must impose 
the identification of points described in (6.15). We will implement this condition in 
the quantum version of the model as a restriction on the physical states. 

We have been assuming so far that the parametrization used for the classical so- 
lutions corresponds to a good set of coordinates in the phase space. In fact, employing 
parallel arguments to those presented in Sec. V for the type II Bianchi model it is 
possible to show that the transformation from the triad and the spin connection to 
the chosen set of parameters is analytic in the considered physical solutions to all but 
the scalar constraint (except at /3 = 0, tt). We will now briefly discuss the main lines 
of this proof. 

Let (gij, Kij) be the metric and extrinsic curvature for the type I Bianchi model, 
which are analytic in the triad and the spin connection for gu a non-degenerate metric. 
We first introduce the matrix Mj = {gi)^M^j, where M^j e SL{3, H) is the matrix 
which appears in (3.2) and leads to a diagonal form for the metric gij, and gi are the 
components of the corresponding diagonal metric. The matrix M satisfies the relations 

detM > , ((M-1)*)/(7pq(M-1)Q, = 5ij , (6.16) 
{{M-y)/KpQ{M-')<^j = XjSjj , (6.17) 

with 

^^ = ^7^-' detM =(^1^2^73)^- (6.18) 
ZI\gi 

One can use Eqs. (6.16,17) to determine M in terms of gij and Kjj. It is convenient 
to adopt the decomposition M = RMt, with R G S0{3), and Mt a lower-triangular 
matrix that includes diagonal elements different from the unity. Eq. (6.16) fixes Mt 
analytically as a function of gjj. On the other hand, the non-diagonal components of 
the system of equations (6.17) provide us with the matrix R as an analytic function 
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of {gij,Kij), as far as A/ ^ Xj for all different I, J = 1,2,3. From expressions 
(6.18) and Eqs. (3.7,8), one can easily check that this is indeed the case if 
{po + 2p_|_ — P-) > 0. Substituting the resulting matrix M, the diagonal components 
of the system (6.17) determine A/ analytically in terms of {gij,Kij). One can then 
identify A/ and detM with their explicit expressions in our parametrization (through 
Eqs. (6.18)), and subsequently compute the parameters pQ,p+,p^ and t as analytic 
functions of {gu, Kjj). Inserting these functions into the expressions (3.7) for gj, 
we obtain the matrix M^j = {gj)~^ Mj, which turns out to depend analytically on 
{gij,Kij). Finally, the parametrization used for j in (6.2) and (6.11) is well- 
defined and analytic, except at /5 = and /? = tt, points that can be interpreted as 
the poles of the two-sphere coordinatized by o; = 2q; and We consider then the 
non-analyticity aX (3 = and = tt as a failure in the coordinatization of S'^. In order 
to go to the reduced phase space, one has only to impose the scalar constraint (3.6) on 
Po, p+ axidp-, and restrict oneself exclusively to positive po (then, {po + 2p+ —p-) > 
is automatically satisfied) . 

VII. Bianchi types I and II: Quantization 

Once we have identified the symplectic structures and the reduced configuration spaces 
for Bianchi I and II, we turn to the task of the canonical quantization of these models. 
Let us start with Bianchi type II. In this case the reduced phase space is the cotangent 
bundle over Cf, , ^ x S^. Keeping the notation introduced in Sec. IV, we parametrize 



by the angle ©, while _|^^ is defined by the constraint Ho = y^+^ + n_^, with 
no,n+,n_ e JR^. a natural set of elementary variables in this reduced configuration 
space is then provided by (11+, n_, c = cosO, s = sinO). As generalized momentum 
variables we choose [9] 




L+ = U+U + -^T, 



(7.1) 







n2 

L- = n_v+ -^T, 



(7.2) 







(7.3) 



where the momenta (T, U, V, Z) are canonically conjugate to (Hq, n_|_, n_, O). Our set 
of reduced phase space variables commute with the constraint (4.15), with respect to 
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the Poisson-brackets structure. On the other hand, the only non-vanishing Poisson- 
brackets among the configuration and momentum variables are: 

{n+,L+} = n+, (7.4) 
{n_,L_} = n_, (7.5) 

{c,Le} = -s, {s,Le} = c. (7.6) 

So, under the Poisson-brackets, the chosen set of variables forms the Lie algebra 
L{T*GL{1,M) X T*GL{1,JR) x E2) where T*GL{1,M) = IR(s)IR+, (s) is the semi- 
direct product, and E2 = 1R^@S'0(2) is the Euclidean group in two dimensions [14]. 
Our next step consists of finding a unitary irreducible representation of the Lie algebra 
(7.4-6). We choose as our representation space the space of distributions '0(n+, n_, ©) 
over _|^^ X S^, and define on it a set of operators (H-i-, n_, c, s, L+, L_, Lg) such 
that their only non- vanishing commutators correspond to i^ times the Poisson-brackets 
(7.4-6). The action of these operators on V'(n+,n_,0) can be consistently defined in 
the form 

(n±v)(n+,n_,e) = n±v(n+,n_,e) , (7.7) 

(^a/')(n+,n_,e) = cosGV'(n+,n_,0) , (7.8) 
(^sV')(n+,n_,e) = sinGV'(n+,n_,e) , (7.9) 
(L+V')(n+,n_,e) = -ihu+-^i;{u+,u_,o) , (7.10) 

(L_V')(n+,n_,e) = -inn_^V'(n+,n_,e) , (7.11) 

(loV') (n+, n_, 6) = -m— ^(n+n_, e) . (7.12) 

The inner product in the space of quantum physical states can then be expressed as: 

<^\il^>= I liH , (7.13) 
where Q = £^ x S^, and the measure /j, is given by 

,= ^^de. (7.14) 

The measure fi is simply the product of the measures that correspond to the uni- 
tary irreducible representations of the groups T*GL{1,M) (twice) and E2 [14]. It is 
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straightforward to check that our operators are self-adjoint with respect to the scalar 
product (7.13). Furthermore, there is no nontrivial subspace in our representation 
space which remains invariant under the action of the operators (7.7-12). This proves 
that the constructed representation of the algebra (7.4-6) is unitary and irreducible. 

Let us consider now the type I Bianchi model. In this case, the reduced configura- 
tion space can be identified with £^ x S'^ x (with certain restrictions still to be 
imposed in the space of physical states that come from the identification of points in 
(6.15)). The sphere S'^ is parametrized by the angles a and f3{f3 e [0,7r),Q; G [0, 27r)), 
the circle 5"^ by the angle 6, and £^ ^-^ is defined by the constraint po = + P-^, 

with p_|_,p_,po £ IR^- From now on, we adopt the compact notation p = (p_|_,p_), 
7 = (a, /3, 9). Given the form of the reduced configuration space, it is natural to choose 
the following over-complete set of configuration variables: 

{p+,P-,c,s,ki,ko,k^i) , (7.15) 

where 

c = cos^, s = sin^, (7.16a) 

ki=Yl{d,P), ko = Yl>{d,p), k-i=Y-\d,P), (7.166) 

and Yi^{d,/3) are the spherical harmonics on the two-sphere [25]. All the functions in 
(7.15) are real, except ki and k-i, which satisfy 

{kiy{a,(3) = -k-i{a,(3). (7.17) 

In this case, our generalized momentum variables are 

L+=p+U+^T, L_=p_V+^T, (7.18) 
Po Po 

La=F,, Lf-^^) = e±^"(±iF2-cot/3Fi), Le ^ , (7.19) 

with {T,U,V, Fi, F2, Fs) the momenta canonically conjugate to (po, <5, 6*). 

Following a similar discussion to that presented for Bianchi II, one can check that 
our set of variables forms the Lie algebra L{T*GL{1,M) x T*GL(1,IR) x x E2) 
(with respect to the Poisson-brackets structure), where = lR^(s)SO{3) is the Eu- 
clidean group in three dimensions [14]. We now proceed to find a unitary irreducible 
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representation of the corresponding Lie algebra of Dirac observables on the space of 
distributions 'ip{p, 7) over the reduced configuration space £^ _^^x S'^ xS^. Parallel to 
the situation in type II, the configuration operators that correspond to the variables 
(7.15) act as multiplicative operators on the chosen representation space, while the 
action of the momentum operators (7.18,19) can be defined by 





(P,7) 




iP,l) 




{P, 7) 




(P,7) 




{P, 7) 



d 



•0(^,7) 



d 



dp 

-i^^^(Pr 

The only non-vanishing commutators in this algebra are 



d_ 
da 



i^ip, 7) 



[n+,L+] 



+ 5 



[n_,L_ 



[c, Le] = -ihs, [s, Le] = ihc, 
[La,k±i]^±hk±i, [Lf^^^ykm] = - m{m ± 1) kr. 



[La, L^a.,13) 



(7.20) 

(7.21) 

(7.22) 
(7.23) 
(7.24) 

(7.25) 
(7.26) 
(7.27) 
(7.28) 



where m = 1,0,-1. The inner product in the space of quantum states, without 
imposing the quantum analogue to the identification of points given by (6.15), takes 
on the expression 



<(f>\'i/j> = / 
Jq 

where Q = x x S^, and the measure /j, is 



(7.29) 



dp+ dp- . ^ ,^ ,^ 
fj, = sm pdpdadO . 



(7.30) 



P+ P- 

Note that all the Dirac observables are represented by self-adjoint operators, except 
km and , which satisfy 

(fc^)t = i-irk-m , (7.31) 
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Let us impose now the identification of points (6.15) in the reduced phase space as a 
restriction on the physical states & L'^ (y^t+ +) ^ '^'^ x , //^ . In the representation 
that we have chosen, the required restriction can be stated as 

3, ^) = — a, TT — /?, ^ + tt) . (7.33) 

The operators that correspond to the canonical variables Fi and F2 in the reduced 
phase space are clearly given, in our representation, by Lg, and Lq, respectively. It is 
easy to check, using Eq. (7.33), that 

(Zai^) (P, «, 13,0)^- (Zai^) {P: -a,n-/3,e + TT) , (7.34) 

and similarly 

(l^V) (P, «, P,0) = - (LoiP) {p, TT - ^ + tt) . (7.35) 

On the other hand, and apart from the factor ordering ambiguities that must be irrele- 
vant in the classical limit, the variable F2 in the reduced phase space can be represented 

by the symmetrized operator corresponding to = — i(e~'"L_|_ +cot P La), where e^" 
and cot (3 must be expressed in terms of km, m — 1, 0, —1. It is then possible to show 
that 

(l^V) (p, «, P,0) = - (l^V') (p, -a,7T-P,e + 7r)+ o{n) . (7.36) 

Eqs. (7.33-36) guarantee that, in the classical limit, the points (6.15) are identified in 
the reduced phase space, recalling that d = 2a and, thus, Fi = 

Therefore, the space of quantum physical states for the type I Bianchi model 
is simply the Hilbert subspace of functions in x x S^,fxj that satisfy 

relation (7.33). Using Eqs. (7.29,30) and (7.33), the inner product in this space is 
easily computed to be 

<0|V'> = 2 / ^ / ^ [ sinPdp [ dd [ de^{p,-f)i;{p,-f), (7.37) 
JjR+ P+ Jm+ P- Jo Jo Jo 

so that we can restrict our attention to the sector 9 G [0,7r), and identify the space 
of physical states with x S'^ x [0, tt),//^. The restriction to this Hilbert 
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space of the set of operators previously defined in [C,^^ have then 

a well-defined action in the quantum physical states for Bianchi type I. 

We have thus succeeded in quantizing the full type I and II Bianchi models fol- 
lowing the non-perturbative canonical approach. The physical interpretation of the 
quantum theories so-constructed is completely analogue to that presented in Ref. [9] 
for the diagonal Bianchi models. The results will be discussed elsewhere [26]. 

VIII. Bianchi Type I: An Alternative Quantization 

Wc want to discuss now a difi'erent approach to the quantization of the type I Bianchi 
model, making use of the symmetries that are present in the scalar constraint at the 
classical level [9]. We will show that the quantum theory obtained in this way is 
equivalent to the reduced phase space quantization of this Bianchi model. 

We will restrict our attention to the sector of positive definite metrics. These 
metrics can be represented by real triads in GL(3, H), which can be uniquely written 
in the form 

E'a = {MTyj{Mnyj,R^''{jL) , (8.1) 

with R^^icfL) e SO{2,), ^l{L = 1,2,3) the associated Euler angles (see Eq.(6.11)), 
and Mt and M^, respectively, an upper-triangular and a diagonal matrix. Let us in- 
troduce the following basis of generators for the upper-triangular and diagonal groups: 

iT')/ = S}Si, iT')/ = S}Si, iT^)/ = SjSi, (8.2) 
(T V = 5/ - 35f 5/ , {T'^)/ ^V3{5}5(-6j5i), (T V = 25/ . (8.3) 

Then, the matrices Mt and can be uniquely expressed in terms of the exponentials 
of these generators as 

MT = nf^ie^^^\ Mo = n^^4e^*^\ (8.4) 

We can consider Xi,i = 1,...,6, and 7l,-L = 1,2,3, as a new set of configuration 
variables. Let us designate their canonically conjugate momenta by and p'^ l, and 
perform a canonical transformation from the Ashtekar variables {Aj, E^) to the new 
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set {xi,^LiP^ iP^ l)- Instead of dealing with p'^'^, it is more convenient to use the 
angular momenta Jl of the Euler angles: 

Jl = (Q-')l • (8.5) 
The matrix Ql^ that appears in this equation is defined by means of the relation [27] 

= e mnR°'pQl , (8.6) 



which implicitly employs the fact that dR^uR^'p is an antisymmetric one- form. 

The spin connection A'j can be obtained by integrating the system of differential 
equations 

{Aj ,E b}-—Q-r-d^--dj^iQ )kl^-^5,5, . (8.7) 

The solution to Eq. (8.7) turns out to be a complicated algebraic expression, although 
the calculations leading to it are relatively simple. We will proceed to discuss the con- 
clusions that can be inferred from the result of these computations without displaying 
the explicit form of Ai°' {xi, ^l, ■, J'^ l) ■ 

We first note that, since the S0{3) connection F/" vanishes in Bianchi type I, 
A/" is purely imaginary. As a consequence, all the momenta {p^,J^l) can be re- 
stricted to be real. On the other hand, substituting given by (8.1-4), and 
^/"(a^j, 7i,p*, J'^l) in Eq. (2.9), the Gauss law constraints for Bianchi type I in the 
introduced set of canonical variables can be rewritten 

Ga = ^eabcR^R'^'j^L eij"- « , (8.8) 
from which we conclude that 

J^L fti . (8.9) 

The vector constraints (2.10) are empty for the type I Bianchi model. We are thus left 
only with the scalar constraint corresponding to (2.11): 

2 2 

cS=(^TrAE) -TtJ^Ae) . (8.10) 
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Substituting {Ai°',E^a) as functions of Xi,^L and in (8.10) (with J'^l set equal to 
zero), we arrive at the expression 

Eq. (8.11) can be considered as a quadratic constraint on the cotangent bundle over 
the six dimensional space coordinatized by the variables Xj. We point out that, from 
Eq. (8.1), this space can be identified with the space of positive definite metrics, since 

= e-^-« (Mt)' AMd)'' i^{Mn)\{MTy p . (8.12) 

The quadratic form that appears in (8.11) endows this space with the natural metric 
[22]: 

ds^ = -6{dx6f + G'^dxidxj , (8.13) 

where Xi = Xi,i = 1, 5, and 

= 6{dx4f + Qidx^f + ^(e-^^^+2^^s((ia;3)^+ 
e-^^^^{dxif + e-^'''-^^'''-{dx2 + xsdxif^ . (8.14) 

The coordinate x^ plays then the role of a time, and de is a global time-like Killing 
vector of the metric (8.13). From Eqs. (8.12) and (8.2-4), the five dimensional space 
coordinatized by Xi is simply the space of positive definite metrics of unit determinant. 
Following now a completely similar analysis to that carried out by Henneaux, Pilati 
and Tcitclboim in Ref. [22], it is possible to prove that this five dimensional space, 
provided with the metric (8.14), can be identified with the coset space ^go{3)^ t^^]" 
Furthermore, they showed that, in the metric representation, the quantum states for 
Bianchi type I can be decomposed in the form 

f{g)= [ dX+ [ dX-fi{X) fda !\in(5d(5 f dej{\^)e^,R^^){g) , (8.15) 

JIR JIR JO Jo Jo 

where A = (A+, A_), /i(A) is a specific measure over IR^, (a, /3, 9) = (7^) are the Euler 
angles that parametrize the matrices R{'Jl) G S0{3), and ex,R(-y){g) are generalized 
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"plane waves" which satisfy the operator constraint associated to the classical con- 
straint (8.11). For any matrix -R(7l) £ 50(3), ex,R{j^){g) is defined as the "plane 
wave" of the rotated metric g: exii(^^){g) = ex,i{RgR^), with e\j{g) given in the 
parametrization (8.2-4) by the expression 

e;^ /(^) = e~''^°^'^^^^+^''^++^)^*+(^'^-"''^^^^ . (8.16) 

The "plane waves" (8.16) are eigenfunctions of the momenta operators P4,P5 and 
P6, defined as self-adjoint with respect to the metric (8.13,14) [22]. The respective 
eigenvalues are A+, A_ and Aq. In order to fulfill the operator constraint corresponding 
to (8.11), (Ao, A4-, A_) must satisfy the relation 

-Ao' + A+2 + A_2 = , (8.17) 

which is a direct analogue of the scalar constraint (3.6). In fact, it is possible to check 
(from Eqs. (8.2-4), (8.12), and (3.4)) that the set of parameters (/3o, /3-) employed 
in Sec. Ill to describe the diagonal metrics for the type I model are related to the 
coordinates (a;"^, x^, x^) by means of the transformation 

Po = ^(2a;^ + x^) + Ci , /3+ = ^(x^ + 2a;^) + C2 , P- = -x^ + C3 , (8.18) 
V 3 V 3 

with Ci,C2,C3 some unspecified constants. The canonically conjugated variables to 
the /3's, {po,p^,p-), used throughout our analysis of the Bianchi type I, can then be 
obtained from (Aq, A+, A_) by completing the canonical transformation (8.18): 

po = ^(2Ao + A+), p+ = -^(-Ao + 2A+), p_ = -A_ . (8.19) 

If we want to consider exclusively inequivalent spacetime geometries, and not only 
different positive definite spatial metrics, we have to restrict the range of po,P+ and 
P- to lie in the positive real axis, as we proved in Sec. VI. Changing coordinates 
from (Ao,A+,A_) to {po,p+,p-), imposing that po = y^p+^ +P-^i and restricting 
p+,P- e IR"*", we arrive, from (8.15), at the following decomposition for the quantum 
sates of the type I spacetime geometries: 

F{g)= dp+dp-fl{p) da sinPdp ddF{p,^)ep^R^^^{g) , (8.20) 
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where p = {p+,p-), fl{p) is certain measure over £^ and ep ji(^y^{g) are the gener- 
alized "plane waves" expressed in terms of p+,P- and 7. The wave function F{p,^) 
characterizes uniquely the quantum state F{g) [22], so that, instead of the metric rep- 
resentation, we can select the representation (p, 7) for the quantization of the model. 
The decomposition (8.20) induces the following inner product: 

<F,G>= / dp+dp-jl{p) da sin^dp / deF*{p,-f)G{p,-f) . (8.21) 
Jc+_^ Jo Jo Jo 

It is obvious that the representation (p, 7) coincides with that used in Sec. VII for 
the quantization of Bianchi I. Furthermore, the inner product (8.21) in the space of 
quantum states differs only in the choice of the measure ji(p) from that determined in 
(7.37). As a consequence, we conclude that the two quantum theories here discussed 
for Bianchi type I result to be unitarily equivalent. 

IX. Conclusions 

We have succeeded in completing the canonical quantization of the type I and II 
Bianchi models while keeping the totality of degrees of freedom of these homogeneous 
gravitational minisuperspaces. Our analysis generalizes the works existing so far in 
the literature [9,10], which had concentrated their attention on the diagonal reduction 
of these systems. 

We have first calculated the explicit expressions of the general solution for these 
two Bianchi types, using the fact that the classical evolution problem can always be 
brought to the diagonal form by a change in the set of left-invariant one-forms on the 
leaves of the homogeneous foliation. This is the first time, to our knowledge, that the 
general solution for Bianchi II has been explicitly displayed. 

The classical solutions have been written in both geometrodynamic and Ashtekar 
variables, restricted to the sector of positive definite metrics. We have determined 
the sets of solutions that correspond to different spacetime geometries, eliminating 
the overcounting of physical states. For both types I and II, the parameters used to 
describe the relevant non-degenerate physical solutions have been shown to define a 
good (analytic) set of coordinates in the phase spaces of these models. The presence 
of the non-diagonal degrees of freedom plays an essential role in the proof of this 
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statement. If we had performed a similar analysis in the reduced diagonal models, the 
phase space coordinates associated to the diagonal degrees of freedom could have been 
extended to a wider range of analyticity, because, in this case, one can consistently 
consider as fixed the set of left-invariant one-forms in the homogeneous foliation. The 
physical states of the quantum theories constructed thereafter for the diagonal models 
can be interpreted as dependent on three-geometries with some preferred directions. In 
our description, however, one is forced to treat the directions in the homogeneous slices 
as being interchangeable, as far as this is allowed by the symmetries of the Bianchi 
model. Apart form the maintenance of the extra degrees of freedom, this is the main 
difference that arises in the study of the full Bianchi types I and II with respect to the 
analysis of their diagonal counterparts. The constraints of the systems are identical 
for both kinds of models, since the whole classical evolution can always be put into 
diagonal form. 

Performing a transformation from the Ashtekar variables to the sets of phase space 
coordinates introduced for Bianchi I and II, we have endowed the reduced phase space 
of each of these minisuperspaces with an analytic symplectic structure. For the type II 
Bianchi model, the symplectic form obtained in this way can be identified with that of 
the real cotangent bundle over the reduced configuration space x S^; for type I, 

the symplectic form corresponds to the real cotangent bundle over x 5^ x 5^, 

with an additional identification of points in S'^ x that we have chosen to impose 
at the quantum level. 

We have proceeded to quantize the models by selecting a complete and closed 
*— algebra of Dirac observables in the reduced phase space of each of these systems. 
We have then constructed an explicit unitary irreducible representation of that algebra 
on the space of distributions in the reduced configuration space. In this way, the 
*— relations in the algebra of observables have been straightforwardly implemented in 
the adopted representation as adjointness relations among the quantum operators and 
the inner product in the space of quantum states has been subsequently fixed. The 
physical interpretation of the mathematical framework obtained here is essentially the 
same as that presented in Ref. [9] for the diagonal reduction of the Bianchi models. 

Finally, we have outlined a different approach to the quantization of Bianchi type I, 
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using the existence of a conditional symmetry [29] to span the quantum states in 
terms of generahzed "plane waves" . This alternative approach may be useful for mini- 
superspace models whose algebra of Dirac observables is not sufficiently known, so 
that the canonical quantization program cannot be yet applied to completion in these 
cases. For Bianchi I, the two quantum theories analyzed in this work have been shown 
to be unitarily equivalent. 

We have thus provided the full type I and II Bianchi models with the framework 
needed to address cosmological problems from the quantum point of view. Our study 
also illuminates the role played by the additional non-diagonal degrees of freedom, 
which turn out to be frozen in the Hamiltonian description. Let us finally point 
out that, from our discussion in Sec. VII, the Hilbert subspace of states that are 
separable as a function in times a function in the "non-diagonal" part of the 

reduced configuration space can be immediately associated with quantum physical 
states for the corresponding diagonal case, by simply neglecting the dependence on 
the non-diagonal configuration variables. One can thus analyze the implications of the 
quantum reduction of degrees of freedom for these two kinds of minisuperspace models, 
as toy models in which one can check the possible validity of the minisuperspace 
approximation in the full quantum theory of gravity [30]. 
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